Correlations between 6>i 3 and 823 in a very long baseline neutrino oscillation experiment 
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A very long baseline experiment, with L/E of 1.6 x 10 4 m/MeV, has been proposed as the 
ideal means to measure #13, including its sign, for neutrino energies less than 50 MeV. Higher 
energies require consideration of the earth MSW effect. Approximating a constant density mantle, 
we examine the v e -v^ and v^—v^ oscillation channels at very long baselines with particular focus 
upon the deviation of #23 from maximal mixing and the sign of #13. It is demonstrated that 
measurements of these oscillation channels in this region are sensitive to 613, including its sign, and 
to 023, including its octant, with a significant correlation between the value of #13 and the difference 
of #23 from maximal mixing. 
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INTRODUCTION 

All existing neutrino experiments (except the LSND 
appearance experiment) can be understood within the 
context of three-flavor neutrino oscillations. By relating 
three mass eigenstates to the flavor states via a unitary 
transformation parameterized in the canonical manner 
0, one has six free parameters with which to fit the 
experimental data: two mass-squared differences, three 
mixing angles, and one Dirac CP phase. Solar neutrino 
experiments [j| along with KamLAND p| point to the 
large mixing angle MSW solution with a mass-squared 



difference A 2 i 



1(T 5 eV 2 (with A 



and mixing angle 6*12 ~ 0.58. The results of the Super- 
Kamiokande atmospheric neutrino experiment |5j , in ad- 
dition to the K2K experiment |(|, indicate another mass 
scale with A 31 ~ 2 x lCP 3 eV 2 and a near maximal mix- 
ing angle # 2 3 ~ 0.79. To date, of the reactor experiments, 
data from CHOOZ provide the most stringent limit on 
the magnitude of the remaining mixing angle | ^13 1 < 0.22 
at 90% CL for A 31 =2x 10~ 3 eV 2 7]. No information 
is known about the level of CP violation in the neutrino 
sector; hence, there are no limits on the phase S. Recent 
analyses can be found in Refs. |a, |9j and recent reviews 

in Refs. OS El 01 El . 

A high-priority goal of the neutrino physics commu- 
nity is to more precisely determine the value of #13. Fu- 
ture reactor experiments have been proposed to improve 
upon the CHOOZ results with the hope of establishing 
a lower bound on the magnitude of this mixing angle 
|l4|. Measuring electron neutrino disappearance from a 
reactor yields a clean measurement of sin 2 #13; however, 
this measurement cannont contribute to determining the 
sign of the angle. In order to reduce the systematic error, 
these experiments will have a near and far detector. An 
increase in sensitivity to sin 2 #13 can be achieved by an 



innovative use of the near detector [l5j . In addition, for 
atmospheric neutrinos traveling through the earth, the 
ratio of multi-GeV e-like to /i-like fluxes is sensitive to 
sin 2 #13 as well as the mass hierarchy [l^. Multi-GeV 
atmospheric neutrinos can also provide information on 
sin 2 #13 and the hierarchy question through a measure- 
ment of the difference between /i + and \x~~ events [l7lll^ |. 
Such a measurement would require an iron calorimeter 
detector such as the MINOS detector [lsj . 

We are particularly interested in measurements that 
are sensitive to terms in the oscillation probabilities 
which are linear in sin #13. It is possible that experiments 
designed to focus upon such measurements might further 
reduce the errors on the value of #13. In addition, other 
information about the mixing angle can be garnered from 
these measurements. We adopt the choice of bounds on 
the mixing angles and CP phase as proposed in Ref . [2(| . 
In contrast to the standard convention, this choice re- 
stricts the CP phase to < 6 < it while extending one 
mixing angle to negative values — ir/2 < #13 < ir/2. The 
advantage is that the parameter space is a connected re- 
gion whenever one takes CP to be conserved. Given this 
choice, it is clear that terms in the oscillation probability 
which arc linear in sin #13 can aid in determining the sign 
of this mixing angle. 

The exact expressions for the oscillation probabili- 
ties necessarily exhibit terms linear in sin #13 |10|: how- 
ever, given that oscillations occur at two different mass 
scales, most current experiments can be understood to a 
good approximation within the context of modified two- 
neutrino oscillations. In Ref. |2l| . one finds perturbative 
expansions about the small quantities #13 and/or the ra- 
tio of the mass-squared differences a := A21/A31 ~ 0.03. 
From these expansions, we note that terms linear in 813 
are supressed by a factor of a. Clearly, this small mix- 
ing angle will be of little consequence whenever this per- 
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turbative expansion is valid. Whenever oscillations due 
to the small mass-squared difference interferes with os- 
cillations due to the large mass squared difference, the 
expansion is no longer a good approximation at low or- 
ders [2l| • In terms of the baseline L and neutrino energy 
E, these subdominant oscillations need to be considered 
whenever L/E > 10 4 m/MeV. 

The sub-GeV sample of atmospheric neutrinos in the 
Super-K experiment contains such baselines || . The im- 
plications of these subdominant oscillations for the anal- 
ysis of atmospheric neutrinos has been thoroughly inves- 
tigated in Ref. j^j. Terms linear in #13, in part, could 
explain a possible excess of e-like events in the Super-K 
data. Earlier work examined this excess but in the ab- 
sence of the linear terms j^. Additionally, in Ref. [24| . 
a full three neutrino analysis of the atmospheric data is 
presented. In this work, the authors assume CP conser- 
vation; as they allow only positive 6*13, they must consider 
two CP phases 5 = 0, it. Relating this to our choice, neg- 
ative #13 corresponds to S = it. Small differences between 
the two choices of the CP phase are a direct result of the 
terms linear in sin #13. Finally, in the analyses presented 
in jH, Ei| , the authors note a small asymmetry in \ 2 
about #13 = attributable to such linear terms. Though 
statistically insignificant, two local minima in \ 2 were 
found, one for positive #13 and one negative. 

In Ref. [25l |. it was determined, in vacuo, that the v^-v^ 
and v^—Vp oscillation channels exhibited signficant linear 
dependence upon #13 whenever one had a baseline-energy 
ratio of L/E = 2(2n+l)vr/Ai 2 . With n = 0, the first A 2i 
oscillation trough, a twenty-five percent effect was found 
on the v^-v^ channel for the presently allowed range of 
6»i 3 . This ideal baseline is L/E = 1.6 x 10 4 m/MeV. 

A purpose of this work is to extend the results of 
Ref. Hi] to include the MSW effect [27| so that one may 
employ higher energy neutrinos to probe the region of in- 
terest cited above. We show herein that even with higher 
energy neutrinos an experiment which is sensitive to lin- 
ear terms in #13 is still beyond present technology but 
perhaps more feasible. We also examine the correlation 
between the value of #13 and the deviation of #23 from 
maximal mixing. We take as known all parameters but 
these two mixing angles; granted, a full consideration 
of all parameters can have significant impact upon such 
correlations j2^. We first consider vacuum oscillations 
which showcase the essential feautres. In the subsequent 
section, we permit matter interactions via the MSW ef- 
fect. An approximate analytical expression is derived to 
demonstrate the relation between the two mixing angles. 
The general features found in the vacuum analysis re- 
main. We then examine through numerical means the 
implications of the results for designing the ideal experi- 
ment that would be most sensitive to terms linear in #13. 
Finally, we note that two existing three-neutrino analy- 
ses [23, might demonstrate the correlation between 
#13 and #23 • 



VACUUM OSCILLATIONS 

We shall assume CP is conserved. This assumption is 
indeed significant as all terms in the oscillation proba- 
bilities which are linear in sin #13 will be modulated by 
either cos 8 or sin<5. Thus if CP is violated, our results 
will be modified at the quantitative level. However, in 
the absence of any evidence for CP violation in the neu- 
trino sector, we will default to the simplest case of no 
violation. 

In vacuo, the relativistic limit of the neutrino evolution 
equation in the flavor basis is 

i9t"f = ^UMrfv f (1) 

for neutrino energy _E, mixing matrix t/, and mass- 
squared matrix M. = diag(0, A21, A31). This first order 
differential equation is easily solved to give explicit ex- 
pressions for the oscillation probabilities. The probability 
that a neutrino of flavor a will be detected a distance L 
from the source as a neutrino of flavor (3 is 

3 

V a() {L/E) = 5 a0 - 4 (UajUakUpkUpj) sin 2 Vjh (2) 

k<j 
j,k = l 

with (fjk := AjkL/4E. As we are interested in the oscil- 
latory region for the small mass-squared difference A21, 
we will assume that the oscillations due to the two larger 
mass-squared differences are incoherent; that is, we take 

(sin 2 (p 3 i) = (sin 2 ip 32 ) = ]r. (3) 

The v e -v e oscillation channel exhibits no linear depen- 
dence on sin #13, so we neglect it in our analysis. Using 
the shorthand notation Sjk := sin(9jfc and Cjk ■= cos8jk, 
the interesting oscillation channels in the prescribed limit 
become 

7% = [ 4 rsin2#i2COs2#i2sin2#i 3 ci3sin2#23 

1 • 2 ca 2/2 2 2 \l • 2 

+ sm 2#i 2 c 13 (c 23 - s 13 s 23 )Jsm cp 12 

+ isin 2 2#i 3S 2 3 , (4) 

= 1 - [cos 2 2#i 2 s 2 3 sin 2 26<23 

+2sin2#i2COs2#i2Si3 sin2#23(c 2 3 — s 2 3 s 2 3) 

1 • 2 nn I 2 2 2 \2i • 2 

+ sm 2#i 2 (c 23 — S13S23) \ sm ^12 

— 2Ci 3 S 2 3(l — c 13 s 23)' (5) 

TV = [- sin 2 26»i 2 (c^ 3 - s 2 3 s 2 3 )(s 2 3 - s\ z c\ z ) 

+ sin 2#i2 cos 2#i 2 si3(l + s 2 3 ) sin 26> 2 3 cos 2# 2 3 

+ cos 2 2#i 2 s 2 3 sin 2 2# 2 3] sin 2 ipi2 

+ ic 4 3 sin 2 2#23- (6) 
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As stated, for the mass-squared difference A31 = 2 x 
f CP 3 , the results of CHOOZ 7j restrict the magnitude of 
#13 by 



1 013 1 < 0.22 (90% C.L.). 



(7) 



Additionally, an analysis [29j which aims to determine 
the octant of 023 indicates a range 

0.67 < 6> 23 < 0.91 (90% C.L.). (8) 

Permitting negative values of #13, the allowed deviation 
of 0i3 from zero is on the same order of magnitude as 
the allowed deviation of 023 from maximal mixing. We 
introduce the parameter e to indicate the octant of 023 
via 



723 



7T 

4 +£ 



(9) 



so that negative (positive) e indicates that 023 is in the 
first (second) octant. 

We expand the in vacuo oscillation probabilities in 
Eqs. JUEIj! to first order in 0i3 and e. We utilize the 
relations 

s 2 3 = -(1 + sin2e), sin 20 23 = cos 2e, (10) 



-23 



= -(l-sin2e), cos 20 23 = - sin 2e. (11) 



Beginning with the v^-v T oscillation channel, we find 
that the term linear in 0i 3 is effectively screened due to 
the near maximal mixing of 023 . The result is a modified 
two-neutrino oscillation probability 



1 



■sin 20i2 sin ¥>2i- 



(12) 



The Super-K atmospheric experiment indicates that v^- 
v T occur most prominently at the A31 mass scale. With 
this in mind, we anticipated this result. 

From Ref. j2{|, one finds that in the oscillation chan- 
nels V efi and terms linear in S13 are not overly 
screened. It is these same terms which are responsible, 
in part, for the possible excess of e-like events in the sub- 
GeV sample of the Super-K atmospheric data (2^ • Drop- 
ping higher order terms, we find that the v^-v^ channels 
exhibit the behavior 



sin"" 20i 



- -e + cot(20 12 )0i3 



For v^—Va oscillations, the probability becomes 



1 



sin 2 20i 



cot(20 



12^13 



sin if 2i- 



(14) 

The effect of 0i3 upon these channels is largest whenever 
the dynamical term sin 2 if 21 is maximal. At this func- 
tions first peak, one has a value of L/E = 27r/A2i ~ 



1.6 x 10 4 km/GeV. This is the result found in Ref. [2g. 
There we suggested that an ideal experiment would be 
to avoid the MSW effect by limiting the energy to a very 
low value, say E < 50 MeV. This is far from possible 
with existing technology, so we expand the discussion to 
include matter effects. 



MATTER OSCILLATIONS 

We consider neutrinos which traverse a constant den- 
sity mantle. By doing so, we avoid the issue of parametric 
resonances for neutrinos which go through the core |3l| . 
A constant density mantle is a relatively good approxi- 
mation; however, it limits our baseline to about 4000 km 
|32j . This, in turn, places a limit upon the neutrino en- 
ergy if we are to reach the peak of the A21 oscillations. 
In what follows, we consider only neutrino energies below 
1 GeV and assume a constant average mantle density of 

3.5 g/cm 3 mm. 

For neutrinos traveling through matter, one needs to 
modify the Hamiltonian to include an effective potential 
which operates on the electron flavor exclusively |27j | . In 
the flavor basis, the neutrino evolution equation in matter 
becomes 



id t vf 



1 

2E 



UMU^ + V 



(15) 



where the operator V operates on the electron flavor with 
a magnitude V — \/2GfN £ , where Gf is the Fermi cou- 
pling constant and N e is the electron number density. 
For a mantle density of 3.5 g/cm 3 , this yields a poten- 
tial around V ~ 1.3 x 10 -13 eV. We note that for anti- 
neutrinos, we need to change the algebraic sign of this 
potential. We shall consider only neutrinos below. 

For sub-GeV neutrinos traversing the earth, matter af- 
fects are most easily addressed in the propagation basis 
developed in Ref. [22j • We will briefly review their deriva- 
tion. In the absence of CP violation, the vacuum mixing 
matrix U can be written as the product of three rotations 



t/ = t/l(023)EM013)EM012) 



(16) 



where Uj(0) is a proper rotation by angle about the 
j-th axis. As U\(ff) commutes with V, we may rewrite 
the evolution equation 



sin 2 v?2i. (13) id t v' = 



-3-(7 3 (012)X(73(012) t 



u 2 (e 13 yvu 2 (d 



13 J 



(17) 

with V 1 = U2(0l3) 1 U 1 (023)^/- 

By conjugating the Hamiltonian in this basis via a lo- 
cally defined U 2 (8), this new propagation basis can be 
approximately described by a Hamiltonian H which is 
block diagonal. This correction to 0i3 is given by 

2 sin 20 13 EV 



tan 29 = 



s? 2 A 2 i - 2 cos 26 13 EV" 



(18) 



As A31 dominates the denominator and 2EV j A 3 \ < 
0.13, one effectively has 

„ sin 29 13 EV 



*31 



This correction results in a modified mixing angle 



'13 



713 



(19) 



(20) 



With this additional rotation, we define locally the prop- 
agation basis with v = U 2 {&)^v' and Hamiltonian H 



H = 







A 3l /2E + s\ 3 V 
where the block is given by 



(21) 



H 2 - — U 3 (9 12 ) 



V 





(22) 

As the correction depends upon the local density, 
we must consider its temporal (spatial) derivative in the 
evolution equation. Using the product rule, we note 



d t v = d t [u 2 (eyy + u 2 {6)^d t v'. 



(23) 



Letting a be the generator of the rotation so that U 2 (0) 
e l9a ', we have 



d t U 2 (9)^ = -iaU 2 (9)^d t 



(24) 



Dropping insignificant terms, one may write the evolu- 
tion equation in the propagation basis as 



id t v = (H + ad t 9)v 



(25) 



Given the neutrino energies and baselines under con- 
sideration in our analysis, the term St9 is negligible. For 
our purposes, it is sufficient to consider the evolution 
equation 



idtv = Hu. 



(26) 



The block H 2 in this Hamiltonian can be easily diago- 
nalizcd in closed form with eigenvalues X±. Of dynami- 
cal relevance is the difference in these eigenvalues which 
yields the effective constant density mass-squared differ- 
ence 



A™ = Aai^ cos 2 2M1 - E/E R ) 2 + sin 2 20 12 , (27) 
where we have defined the resonance energy to be 

A 2 i cos26>i2 



E, 



2Vcf 3 



(28) 



Fixing the solar mixing angle 9\ 2 = 0.58, we find the 
resonance energy in the mantle is Er ~ 120 MeV. Mat- 
ter effects require an accommodation to the other mass- 
squared difference A^ = A31 — 2_EA_. The mixing angle 
which achieves this diagonalization satisfies 

sin 20i2 



sin 20^ = 



cos 2 20i 2 (l - E/E R ) 2 + sin 2 20 i2 



(29) 
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FIG. 1: Ideal baseline for mono-energetic neutrinos traveling 
through a constant density mantle. 



With this diagonalization, we see that, in the flavor 
basis, Eq. $Z$ft becomes 



2E 



(30) 



with the modified mass-squared matrix A4' = 
(OjAqijAIji). The new mixing matrix accounts for the 
matter mixing angles 



v = (7 1 (0 23 )c/ 2 (0r 3 )c/3(0r2) 



(31) 



We now discuss these modifications within the context 
of incoherent (coherent) oscillations of the large (small) 
mass-squared difference. As the matter potential V is 
larger than the kinetic term A 2 i/2E for the energies of 
pertinence, the effect on the smaller mass-squared dif- 
ference is most profound. Examining Eq. I|27() . we note 
that for neutrinos of energies below twice the resonance 
energy, the effective mass-squared difference, and thus 
frequency of oscillations, decreases. For neutrino energies 
greater than twice this resonance energy, the frequency 
increases. Assuming a mono-energetic neutrino source, 
the first peak of sin 2 ip^[ occurs at a baseline 



L = 2nE/A™ 1 . 



(32) 



Setting 0i 2 = 0.58, we see in Fig. JQl how this ideal base- 
line changes as a function of energy. Recalling that the 
resonance energy is 120 MeV, we find that at this energy 
a baseline around 2000 km will yield maximal sensitivity 
to the terms linear in 0^. At twice the resonant energy, 
the necessary baseline has increased to 3700 km. Finally, 
the maximal baseline permitted by our constant density 
approximation is 4000 km, which corresponds to a neu- 
trino energy of 270 MeV. 



Regarding the larger mass-squared difference A31 , the 
modification of this quantity needed to account for mat- 
ter was the subtraction of 2E\_ . This term is at most on 
the order of 10~ 4 eV 2 ; therefore, our assumptions regard- 
ing the incoherence of these oscillations is unaffected. In 
short, we may still take 



(33) 



With respect to the mixing angles, the mixing angle 
6*12 undergoes the greatest degree of modification to ac- 
commodate matter effects. This discussion is given in 
the next section. The effect upon #13 is quite small. Es- 
timating the maximum value of the additive correction, 
we have 



2EV 



*31 



013 < 0.09 013. 



(34) 



This represents at most ten percent correction for neutri- 
nos of energies around 1 GeV. Should we only consider 
neutrino energies below 300 MeV, the correction becomes 
less than three percent. We shall consider this correction 
in our numerical analysis but drop the term in the analyt- 
ical discussion. Incidentally, this correction to 6*13 is the 
only place where the sign of A31 has any effect. Hence, 
we see that the hierarchy ambiguity plays a very small 
role in these calculations. Finally, the remaining mixing 
angle 2 3 is insensitive to matter effects. 

The upshot is that Eqs. are still valid for neu- 

trinos traveling through a constant density medium pro- 
vided we make the appropriate substitutions for the mix- 
ing angles and small mass-squared difference. Neglecting 
the small correction 9, the first order expansion of these 
probabilities about #13 and e have not changed in form 



V, 



efj, 



sin 2 26» 



12 



\-e + cot(20™)0 13 



sinV™, (35) 



and for v^—v^ oscillations, the probability becomes 



- - sin 2 29 



12 



1 

A'" 



cot(26T 2 )f? 13 



sin 2 <^i. 



(36) 



DISCUSSION 



Aside from the adjustment in the ideal baseline, the 
change in the 6>^, mixing angle is the most critical mod- 
ification needed to account for matter effects. In the 
approximate formulae in Eqs. (|35l) and (|36[) . the oscilla- 
tory terms are modulated by sin 2 29™ 2 thereby affecting 
the amplitude of oscillations. Recalling Eq. (I29|l . we note 
that, relative to vacuum oscillations, this factor enhances 
the oscillations in matter provided that the neutrino en- 
ergy is less than twice the resonant energy. The new mat- 
ter mixing angle is maximal whenever E ~ Er and equal 




FIG. 2: The solid line is the ratio sin 2 IQfil sin 2 20 i2 . The 
dotted line is the curve cot 2#53. We have set #12 = 0.58 and 
plot the curves as functions of the dimensionless ratio E/Er. 



to the vacuum value whenever E *C En or E ~ 2Er. For 
neutrino energies greater than twice the resonant energy, 
the factor sin 2 28^ decreases significantly, effectively sup- 
pressing the amplitude of oscillations. In Fig. we plot 
the ratio sin 2 29\ n 2 / sin 26>i2 as a function of neutrino en- 
ergy in order to demonstrate the enhancement and sup- 
pression of the oscillation amplitude relative to vacuum 
oscillations. At an energy of 270 MeV, which is 2.25-Er, 
we find that the amplitude will be suppressed by a factor 
of 0.9 relative to vacuum oscillations. 

Another notable feature of the first order expansions 
of the oscillation channels in Eqs. (|35|l and (|3(j|) are the 
terms which appear within the square brackets. Given 
the lack on knowledge in the true values of #23 and 6*13, 
it is apparent from the equations that a measurement at 
these baselines of, say, the Ve—v^ probability does not di- 
rectly resolve the values of these mixing angles. Indeed, 
we see that a host of parameters can satisfy a measure- 
ment provided the linear relation 



e - cot(26>S)0i3 = const 



(37) 



holds. 

As the slope of this relation, cot 26™ 2 , varies with neu- 
trino energy, its value is represented by the dotted curve 
in Fig. (J2J. For energies greater than the En, the slope 
increases linearly with neutrino energy. Of particular im- 
portance are two points. First, at the resonance energy 
we see that both oscillation channels do not depend lin- 
early on 0i3j a measurement of such neutrinos could at 
best improve the bounds upon #23 while potentially re- 
solving the octant of $23. As our constant density ap- 
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proximation is really only valid up to 270 MeV, we note 
from the graph the value of cot 29^ is 0.54 at this energy. 
A measurement at this energy would have implications 
for both the sign of 6* 13 and the octant of #23- 

The aim of this paper is to determine how a measure- 
ment at this baseline of either the v^-v^ or v^—v^ oscilla- 
tion channel will affect the values of #23 and 6*13. To most 
explicitly demonstrate this, we will create some mock ex- 
perimental data and combine it with the present knowl- 
edge of these mixing angles. To this end, we introduce a 
X 2 function made up of three parts 



2 2,2,2 

X =X13+X23 + XP- 



(38) 



For the first term, we include information from CHOOZ 
via the simple relation 



xi 3 



?13 



9? 3 ) 2 



(39) 



'13 



where the best fit value is 9® 3 = and the standard 
deviation is (J13 = 0.15. Similarly, we include information 
on #23 m the same manner where the best fit value is 
maximal #23 — 0.785 and the standard deviation is 023 — 
0.07. 

The final term in \ 2 represents our mock data. For 
concreteness, we choose the electron-muon oscillation 
channel. As with the others, we define Xt> with respect 
to a best fit P° and a standard deviation which we take 
to be ten percent of the best fit 



xl 



13 7 "23 ) 



(40) 



We assume, as well, a mono-energetic source of 270 MeV 
measured at the ideal baseline of 4000 km. The question 
remains as to what value we assign our measured best fit . 
At present, the only experimentally observed neutrinos 
of this energy and baseline occur in the sub-GeV sample 
of the Super-K atmospheric experiment 0|. A detailed 
analysis of the oscillation channels of interest has been 
performed in Ref. \22\ . Our first order approximations in 
e and #13 exhibit some of the significant behavior given in 
this reference. Briefly, assuming neutrino oscillations and 
CP invariance, one would expect the detector-to-sourcc 
ratio of atmospheric electron-like neutrinos to be 



K r , = Ve. 



rV, 



(41) 



where r is the ratio of muon to electron neutrinos in the 
upper atmosphere. For energies of a few hundred MeV, 
one has r ~ 2. For those which traverse the mantle only 
with baselines over 3 x 10 3 km, we may use the expression 
from Eq. I|35|l along with a first order expansion of V ee 
to determine the approximate expectation ratio 

TZ e ~ 1 + 2 sin 2 29™ 2 [cot(26»™ )6>i 3 - e] sin 2 tp™ . (42) 



In-depth analyses of Super-K |22j, |23J, |29j cite the os- 
cillations discussed herein as a potential explanation for 
an excess of electron-like events in the sub-GeV sample. 
If this excess is indeed real, then our approximation in- 
dicates 



cot(20S)0i3 - e > 0. 



(43) 



Though this expression is not the result of a detailed anal- 
ysis of the experiment, its qualitative features should still 
be valid. With that said, the inequality in (|4^|l indicates 
that the potential excess of electrons in the sub-GeV 
sample can be accounted for by a range of parameters. 
Should #13 = for instance, we find that e < which 
indicates that #23 is in the first octant. Alternatively, 
should 2-3 mixing be maximal, the above inequality in- 
dicates that #i3 is positive. 

This qualitative understanding is confirmed in more 
detailed analyses of experiments which have very long 
baselines. In Ref. |29j . preliminary evidence indicates 
that #23 lies in the first octant. Within this full three- 
neutrino analysis (labeled A), the best fit mixing angles 
are found to be 



'12 



0.575, 



'13 



0.0, 



9 A 
7 23 



0.745, 



(44) 



which in the present context is expressed by e = —0.04. 
Additionally, in Ref. [3o| . a full three- neutrino analysis 
of the world's data is performed. When 2-3 mixing is 
assumed maximal, the best fit mixing angles (labeled B) 
are determined to be 



1 B 
'12 



0.591, 



1 B 

'13 



0.095, 



0g = 0.785. (45) 



These two parameter sets are consistent with the corre- 
lation between #13 and e which we have indicated above. 
Employing the average of the two values for the solar 
mixing angle #12 = 0.58, we find 



7 23 



7 2:S 



'13 



iA 
'13 



= cot 26»i2, 



(46) 



to four percent. Thus, although a small effect, the con- 
tribution to the e-like sub-GeV excess for atmospheric 
data [23 introduces a correlation between the extracted 
values of #13 and e. 

In Ref. [29j it is shown that future atmospheric data 
will further constrain the value of #23 and perhaps de- 
termine its octant, this for #13 = 0. We have here im- 
plicitly generalized this result. Future atmospheric data 
will provide information on the value of the function 
cot(26 ) "2) #13 — e in addition to the cleaner measurements 
of e. 

Though these values are not statistically robust, we 
may use them as an indication of where a best fit point 
might lie. For the best fit, we choose the average value 
for the solar mixing angle and the values for #13 and #23 
from parameter set A in Eq. (|44|l . We may now determine 
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FIG. 3: Given the mock data for the V RV . measurement at the 
ideal baseline for mono-energetic neutrinos of 270 MeV, we 
construct a hypothetical allowed region for #13 and #23- We 
draw the 68%, 90%, and 95% CL contours. 

the acceptable regions for these two mixing angles given 
our mock data and the existing knowledge. From the x 2 
function, we are able to plot the 68%, 90%, and 95% CL 
allowed regions in Fig. ©. The best fit value is 9i 3 — 
0.04 and 6*23 = 0.775. The slope of the upper portion 
of the closed contours in the figure is given by cot 29^ . 
We note that this mock data indicates a preference for 
positive 6*i3 and first-octant 623; however, the possibility 
remains that they lie elsewhere. 



CONCLUSION 

We have expanded the oscillation probabilities V ep , and 
P w to first order in #13 and e, where e indicates the de- 
viation of the mixing angle 623 from maximal. This was 
done for neutrinos propagating in the vacuum and in a 
constant density mantle. For the oscillations in matter, 
it was shown that the amplitude of the oscillations is en- 
hanced relative to vacuum for neutrino energies less than 
twice the resonant energy; the amplitude was suppressed 
at energies greater than this mark. Terms linear in #13 
were of significance whenever the matter-modified solar 
mass-squared difference A 2 ^ was in the oscillatory region. 
The ideal baseline, defined by the first peak or trough of 
these oscillations, was determined as a function of neu- 
trino energy. For a maximum baseline of 4000 km, this 
energy was found to be 270 MeV. At this energy and ideal 
baseline, the relation between the deviation of #23 from 
maximal mixing and #13 was examined by creating mock 
data. The mock data demonstrated what was found in 
the approximate equations; that is, #13 and e are linearly 



related via a slope of cot 26*™2- The energy dependence of 
this slope was elucidated. For an energy of 270 MeV, it 
was shown to be 0.54; this factor increases linearly with 
neutrino energy, assuming the energy is greater than the 
resonant energy. 

It is seen from the mock data that this ideal measure- 
ment could help determine the algebraic sign of #13 and 
perhaps the octant of #23- For a measurement most sen- 
sitive to the octant of #23, one would like the neutrino 
energy to be near the MSW resonance of 120 MeV. At 
this point, the slope is minimal because cot 26*^ = 0. 
To best determine the sign of #13, one would like the 
slope relating the mixing angle and e to be as large as 
possible without completely suppressing the oscillation 
amplitude. To achieve this, one would require higher en- 
ergy neutrinos and, therefore, longer baselines. Our con- 
stant density treatment would no longer be valid as these 
neutrinos would certainly traverse the core of the earth 
and parametric resonances would be of consequence. We 
leave a detailed analysis of this case to future work. 
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